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We show that the RenormaHzation Group formahsm allows to compute with 
accuracy the zero temperature correlation functions and particle densities of 



t;;^- ' quantum systems. 
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O ' In [0] we showed that the Renormalization Group is an efficient tool to compute the 

^ ' ground state and first excited energy levels of a one quantum particle system. The method 
is the following: instead of computing directly the effective potential Vo(xo) defined as a 
^ ' path integral on the Fourier modes: 
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-p (-/?^(xo)) = / n exp (-1^.) , (1) 

' 

I we compute recursively a running potential by integrating each modes after the other. Note 

that the Feynman path integral is considered here with a discretized time t„ = = ne, 

with an arbitrary large number, and n = 0, . . . , + 1. The Fourier decomposition of a 

periodic path contains only a finite number of Fourier modes 
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X{tn) = Xo + -j===Y,^W{i^mtn)Xm + CO., (2) 
V^v + i 

where J2' is from 1 to ^ if is even and from 1 to if is odd. The Xm are the Fourier 
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modes and co'i. = — The discrete action is: 
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SM = eY.Mui\x^\' + eY: VMtr,)) (3) 
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In the Local Potential Approximation |T|], p[, the potential at scale m-1 is given by the 
Wegner-Houghton equation: 



Vm-l{xo) = Vmixo) + ^ log(l + ^S^^) + O(^) (4) 



In the present work we consider only the zero temperature case and neglect as a consequence 
terms of order Our experience tells us that the finite temperature case is more difficult 
to handle with this equation |Q and gives less good results than the Feynman-Kleinert 
variational method In order to use this method to get the correlation functions we have 
to compute the partition function in the presence of a source term: 

. - / ru(f § +v(§)-i(u)§(u)) 

Z{3)= V^] ' (5) 



Note that the expression is written in the continuum only for convenience. As in we 
define recursively the effective potential V^_i{xq) (which is now j dependent) at scale m — 1 

for a constant background path xq = ^ J dtx (t) , as a path integration on the modes Xm and 

^ 

X—m- 

At the end of the ffow of V^{xo), formula (|^) reduces to a classical partition function with 
a source term: 

f dx^ p(v^(xo)-joxo) 

J 2ti(5 

The gaussian integration in (|]) leads to: 

or in terms of the Fourier transform of the source: 

Vi-,ixo) = Viixo) + ^ log(l + }^^^) - ^-^-g (8) 
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By neglecting the source term dependence in the right hand side of equation ( 0), we obtain 
the following approximate relation: 

where now V^(xo) is the running potential without source term defined in equation (^). It 
will be shown below, when computing the particle density, that this approximation is quite 
good. Note that the double sum in the left hand side of (|) is of order so that the whole 
source contribution is of order 

After iteration, (^) leads to the effective potential Vo'(^o) in terms of the effective poten- 
tial without source Vo(a;o): 
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P m=l UJ^ + Vrn [Xq) 

where for convenience we have replaced the double sum by integrals. It's then easy to 
compute the correlation functions using the following relation: 

Since we consider only the zero temperature limit, the integral in (|ll|) is dominated by the 
minimum of the effective potential V [xq) so that: 

{x{h)...x{Q) = -^ -^e'r.^^ (12) 
As an example the two point correlation function is easily deduced: 
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]_ 1 2 ^iWm{t-t') 



In a similar manner we will now compute the particle density which is defined as: 

p{xa) = f dx (x |e-^^| x) 5{x - X,) = ^ \Mxa)\^ e-^^", (14) 



In the limit /3 — > cxo: 



(15) 



where ipoi^a) is the ground state wave function. In terms of the path integral the density is: 



I?§(U)(5(§(r) -§H)exp 



-§ +V(§(U)) 



(16) 



valuable for any value of r. These expression can also be written: 

dxQ f dk 
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Setting j (t) = —ikS{t — r) allows to use the previous formula (|TD|) to directly write: 



dXn f dk f,2^2/2+ikXa^-l3Voixo) 
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where: 



13 ^iul+Vr{xo) 
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after integrating out k, we obtain: 



dxo e" 
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So that in the limit /5 — »• oo: 

p{x 



f ^27ra2(xo) 



-/3Vb{xo) 
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-(a:-xo)V2a^(^o) 



(21) 

^27ra^(xo) 

is a gaussian distribution. The procedure is now very simple. We just have to determine 
the flow of the second derivative of the effective potential. Note that the distribution (^) is 

JV 

the same as the variational Kleinert's one except that for him a^(xo) = J2m=o muj^ +n{xo) 
where Q is the variational parameter 0]. 

As an example consider the anharmonic oscillator. In [|I|] we showed that the running 
potential at scale m is best approximated by a polynomial of order six: 



VM =9l+ '-fx' + 'fx^ + 'fx' 



(22) 



where the initial potential is chosen to be: 



N 1 2 240 4 
Vn{x) = -Xo + — Xo 



(23) 



Since the minimum of the effective potential is located at = it is then easy to see that 



a2(0) 
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. From the properties of the correlation function we can also deduce 



that the first excited energy level is given by: Ei — Eq 
From 0] we have: 
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1.4982 


1.5313 


1.50497 


5.3368 


5.3482 


5.3216 


0.1472 


0.125 



A comparison between the RG density and the variational density is shown in figl. The 
RG density is very near the exact density given in |^ . In particular the RG computation 
greatly improves the value of the density near a; = 0, (p_rg(0) ~ Pez(O)). 

The double well potential case is much more involved. It is well known that the effective 
potential at T = must be a convex quantity. For a very deep potential the truncation of 
the running potential is no more a good approximation of the true effective potential near 
zero, which is almost fiat. Indeed the truncation of the running potential does not allow 
a positive term in the logarithm of the Wegner-Houghton equation and the approximation 
breaks down. We have considered the case A = 2.4 where an expansion at the order 10 is 
necessary to get valuable results: 
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2.4 


0.46498 


0.549 


0.4709 


0.8285 


1.035 


0.7677 


3.48 


1.03 



An another problem with our approach is that a convex effective potential is incompatible 
with the two picks of the true particle density |0] (it's clear that the particle will spend 
most of his time around the two minima of the potential energy) . In statistical mechanics 
language each spin fluctuates most of the time around two values (like Ising spin) . In fact our 
approximation which relates the density to the effective potential breaks down in this case. 
Nevertheless, the behavior of our RG density is quite correct near zero (pkg(O) ~ PexiS^)) 
and for large values of x. Moreover it is much more better than the variational computation 



(fig2.). To get a better result we must go beyond the LPA approximation, and make the 
computation for an arbitrary background path configuration. 
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FIGURES 
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FIG. 1. The RG density of the anharmonic osciUator potential as compared to the variational 
density (dashed line). 




FIG. 2. Comparison of the particle density of the double well potential between the RG den- 
sity and the variational density (dashed line). The RG density is globaly much better than the 
variational one (to compare with exact result see [Q]) but our approximation is not good enough 
to show the double picks structure of the true density (see text). 
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